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Abstract 

Investigations on the internal one-particle density matrix in the case of Bose-Einstein condensates 
with a finite number ( ) of particles in a harmonic potential are performed. We solve the 
eigenvalue problem of the Pcthick-Pitacvskii-typc internal density matrix and find a fragmented 
condensate. On the contrary the condensate Jacobi-type internal density matrix gives complete 
condensation into a single state. The internal one-particle density matrix is, therefore, shown to 
be different in general for different choices of the internal coordinate system. We propose two 
physically motivated criteria for the choice of the adequate coordinate systems which give us a 
unique answer for the internal one-particle density matrix. One criterion is that in the infinite 
particle number limit { N = oo ) the internal one-particle density matrix should have the same 
eigenvalues and eigenfunctions as those of the corresponding ideal Bose-Einstein condensate in the 
laboratory frame. The other criterion is that the coordinate of the internal one-particle density 
matrix should be orthogonal to the remaining {N — 2) internal coordinates, though the {N — 2) 
coordinates, in general, do not need to be mutually orthogonal. This second criterion is shown to 
imply the first criterion. It is shown that the internal Jacobi coordinate system satisfies these two 
criteria while the internal coordinate system adopted by Pethick and Pitaevskii for the construction 
of the internal one-particle density matrix does not. It is demonstrated that these two criteria 
uniquely determine the internal one-particle density matrix which is identical to that calculated 
with the Jacobi coordinates. The relevance of this work concerning a-particle condensates in nuclei, 
as well as bosonic atoms in traps, is pointed out. 
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I. INTRODUCTION 



The second 0+ state of ^^C which is located near the 3a breakup threshold and is called 
the Hoyle state [1] is well known as one of the mysterious 0+ states in light nuclei. The 
understanding of its structure has been actually one of the most difficult and challenging 
problems of nuclear structure. Its small excitation energy of 7.66 MeV is clearly not easy 
to explain by the shell model and, in fact, even the most advanced modern shell model 
approach known as the no-core shell model, fails by far to reproduce its position [2]. 

About 40 years ago Morinaga proposed to assign the 3a hnear-chain structure to this 
Hoyle state [3] . However the observed reduced a decay width of this state which is larger 
than the Wigner limit value was pointed out to be contradictory to the linear-chain struc- 
ture which can give at most only one third of the Wigner hmit value [4] . The large a-decay 
reduced width of the Hoyle state in the ^Be{Of)+a channel was reproduced by a full three- 
body calculation with a semi-microscopic a-a interaction [5], namely an OCM ( orthogo- 
nality condition model ) [6] calculation for the 3a system. This 3a calculation contradicted 
the 3a chain structure of the Hoyle state and gave as the dominant component of the Hoyle 
state the ^Be(Oi') -I- a structure with relative S wave between ^Be(Oi") and a. Since ^Be(Oi") 
consists of two a clusters weakly coupled in relative 5" wave, the Hoyle state was concluded 
to have a weakly coupled 3a structure in relative S waves with large spatial extent, which 
is a gas-like structure of a clusters. A few years later, this understanding of the structure 
of the Hoyle state was reported to be reproduced by full microscopic 3a calculations by two 
groups, namely Kamimura and his collaborators [7] and Uegaki and his collaborators [8]. 
These calculations nicely reproduced not only the position of the Hoyle state but also other 
experimental properties including inelastic electron form factor and EO and E2 transition 
properties. Other states of ^^C below and around the Hoyle state were also well described. 
The microscopic 3a model treatments of ^^C of Refs. [7, 8] were later extended to reaction 
theory in Ref. [9]. Furthermore investigation by antisymmetrized molecular dynamics [10] 
and that by fermionic molecular dynamics [11] starting from a realistic bare N-N force, 
both of which do not assume alpha clustering, have now reproduced all the salient features 
of the Hoyle state [1]. Let us also mention that for a product state of only a few bosons, 
the phase fluctuates. However, we use the word 'condensate' in the same 'loose' sense as 
nuclear physicists are used to say that 'nuclei are superfluid', in spite of the fact that there 
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is not a macroscopic number of Cooper pairs in nucleus. 

Almost 30 years after the first proposal of the ^Be(Oi") + a. structure for the Hoyle state, 
this state was reconsidered in a new light in Ref. [12] involving four members of the present 
authors. They presented the following new model wave function $^(30;), called the THSR 
wave function: 

$B(3a) = ^{exp[- + + X^)] <^{a^)<^{a^)4>{a^)} (1) 
= exp(--|j^^)^{exp(-^^? - 0(ai)0(«2)0(a3)}, (2) 

= Xi - \{X^ + Xs), ^2 = ^2 - X3, = ^(Xi + X2 + X3), (3) 

where (^[ptj] and Xj stand for the internal wave function and the center-of-mass coordinate 
of the j-th a. cluster, respectively, and A is the antisymmetrizer of the nucleons. As shown in 
Eq. (2), the THSR wave function can be regarded as expressing the cluster structure where a 
^Be(Oi")-like cluster ^{exp(— (1/S^)^2)0('^2)0(q;3)} and the cluster couple via an S-w&ve 
inter-cluster wave function exp(— (4/3-B^)^i). On the other hand, Equation (1) shows that 
the THSR wave function represents the state where three a clusters occupy the same single 
O^-orbit exp(— (2/i?^)X^), namely a 3a condensate state which is the finite size counterpart 
of the macroscopic ct-particle condensation in infinite nuclear matter at low density [13]. 
What the authors of Ref. [12] proposed was that the ^Be(Oi") + a structure of the Hoyle 
state can be regarded as being a 3a condensate state and that one can expect in general 
the existence of na condensate states in the vicinity of the na breakup threshold in 4n self- 
conjugate nuclei [14]. It was soon discovered [15] that the microscopic 3a wave functions of 
both of Refs. [7, 8] have overlaps of more than 95% with a single THSR wave function with 
a large size parameter B, implying small overlaps between individual a clusters. 

This above-mentioned striking fact reported in Ref. [15] means without doubt that the 
Hoyle state structure has a strong relation with the a condensation physics in dilute infinite 
nuclear matter. One of the important tasks for the clarification of this relation is to study 
the magnitude of the component of the ideal Bosc-Eiustciu condensation of structureless 
a particles which is contained in the Hoyle state wave function. This was investigated by 
the authors of Refs. [16, 17] who solved the eigenvalue problem of the internal one-particle 
density matrix of the Hoyle state wave functions where the center-of-mass coordinate is 
eliminated. As the internal coordinates for calculating the internal one-particle density 
matrix p they used the Jacobi coordinates, namely and ^2 Eq. (3). They obtained a 
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maximum eigenvalue larger than 0.7 for the normalized density matrix p {Tr p — 1 ), which 
means that the corresponding eigenfunction which is a 05'-orbit is occupied to more than 
70% by the three a particles. This is a large percentage corroborating the almost ideal a- 
particle condensation nature of the Hoyle state, that is the Hoyle state is describablc to very 
good approximation by a product state of these bosons, each in the same 05* orbit [16, 17]. 

For the discussion of the Bose-Einstein condensation of a confined macroscopic system, 
one uses in general the one-particle density matrix in the laboratory frame where the coordi- 
nate system consisting of individual particle coordinates is adopted. However, Pethick and 
Pitaevskii [18] (PP) proposed to use the internal one-particle density matrix by eliminat- 
ing the ccntcr-of-mass coordinate when one discusses the condensate fraction of the system 
where only the centcr-of-mass degree of freedom is excited but all relative degrees of freedom 
are kept unchanged. Their proposal was to point out that the conclusion of the paper by 
Wilkin et al. [19] does not appropriately reflect the physics. Wilkin et al. discussed the 
lowest excitation of a condensate of cold bosonic atoms with attractive interactions which 
rotates with its center of mass in a harmonic potential but keeps all the internal degrees 
of freedom as in its ground state which shows complete condensation. It was found that 
the corresponding one-body density matrix in the laboratory frame has many eigenvalues 
of comparable size, and thus the system should be characterized as a fragmented conden- 
sate [19]. On the other hand PP claimed that if one uses an internal one-particle density 
matrix by eliminating the center-of-mass coordinate it should give a single eigenvalue of the 
order of the number of particles indicating the non-fragmented condensate character of the 
internal part of the system. Unfortunately PP presented only their idea and did not demon- 
strate explicit results of the diagonalization of their definition of the internal one-particle 
density matrix. 

In this paper, following the definition of PP, we construct the internal one-particle den- 
sity matrix of a many-boson system in a harmonic trap and then give the explicit analytical 
form of the eigenvalues and eigenfunctions of this density matrix. We will see that the 
eigenvalues are fragmented, which is contrary to PP's initial objective. We will discuss that 
this result does not mean that the original idea of PP is incorrect but it means that the 
internal coordinate system which PP adopted is inadequate. Actually when we construct 
the internal one-particle density matrix by using the internal Jacobi coordinates, the result- 
ing eigenvalues show complete condensation. We will give two criteria for the choice of the 
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adequate coordinate systems which will give us a unique answer for the definition of the 
internal one-particle density matrix. One criterion is that in the infinite particle number 
limit { N = oo ) the internal one-particle density matrix should have the same eigenvalues 
and eigenfunctions as those of the corresponding ideal Bose-Einstein condensate in the lab- 
oratory frame. This is in line with the general wisdom that in the thermodynamic limit all 
macroscopic quantities shall be the same, regardless whether considered in the laboratory 
frame or in the internal frame [20] . The other criterion is that the coordinate used in the in- 
ternal one-particle density matrix should be orthogonal to (or maximally independent from) 
the remaining {N — 2) internal coordinates, though the {N — 2) coordinates, in general, 
do not need to be mutually orthogonal. This second criterion is shown to imply the first 
criterion. We show that the internal Jacobi coordinate system satisfies these two criteria, 
while the internal coordinate system adopted by PP for the construction of the internal 
one-particle density matrix does not. Furthermore we argue that these two criteria uniquely 
determine the internal one-particle density matrix which is the same as that calculated with 
the Jacobi coordinates. The results of this paper justify the use of the Jacobi coordinates 
in Refs. [16, 17] and, hence, corroborate the almost ideal a-particle condensation nature of 
the Hoyle state. 

The present paper is organized as follows. In Sec. II, we formulate the internal one-particle 
density matrix with respect to the internal wave function of a Bose-Einstein condensate with 
finite particle number in a harmonic trap. Then, the analytical form of the eigenvalues and 
eigenfunctions is presented for the internal density matrix and numerical eigenvalues are 
discussed. In Sec. III. we propose two criteria for the choice of proper coordinate systems for 
internal one-particle density matrix. Finally, the summary is given in Sec. IV. Appendix A 
serves to present the analytical solution of the eigenvalue problem of the density matrix and 
Appendix B is given for the explanation of the bosonic symmetry of the Jacobi-type internal 
one-particle density matrix. The original version of this paper was given in an article on the 
arXiv [21] and a short version of this paper is reported in a letter paper [22] with special 
attention to the cold atom community. 
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II. INTERNAL ONE-PARTICLE DENSITY MATRIX 



First we consider the one-particle density matrix in the laboratory frame for an ideal 
Bose-Einstein condensate with N spinless bosons in a harmonic potential. The result is 
trivial but instructive for studying the nature of the internal one-particle density matrix, as 
will be discussed later. 

The A^-particle Hamiltonian in laboratory frame is presented as 

N 1 N 

^ = E^p? + E^--'^^ (4) 
i=i ^'"^ i=i ^ 

The normalized ground-state wave function of this system is expressed as a product of 
identical Gaussian single-particle wave functions, i.e. 

A /2z/\3/4 



i=i ^ ^ 

where v — muj/2h, and {T*i}^i denotes the set of the coordinates Vi {i — 1- ■ ■ ,N). The one 
particle density matrix in the laboratory frame is defined as 

PiiU^' O = / dr2dr^ ■ ■ ■ drN^*{r, rs, Tg, • • • , ■rjv)$(r', Ts, rg, • • • , r^v), (7) 

= (^^y\^p[-u{r' + r'')]. (8) 

It is noted that the density matrix is independent of the number of particles N and is 
separable with respect to r and r'. The separability originates from the fact that the 
Hamiltonian is separable with respect to different particle operators, Vi and {i — 1, ■ ■ ■ , N), 
in Eq. (4), or equivalently from the fact that the wave function of Eq. (6) is separable with 
respect to different particle coordinates, {i — 1, - • • ,N). 

The nature of the single particle orbits and their occupation probabilities in the relevant 
system can be obtained by solving the eigenvalue problem of the density matrix, 

j ptl{ry)<t>{r')dT'^\4>{r), (9) 

where 0(r) is the single particle orbit, and A is its occupation probability. This equation can 
easily be solved, and we find that the density matrix has only one non zero eigenvalue A = 1 
with one eigenf unction, namely, the zero-node S-w&Ye Gaussian 0(r) = (2i//7r)^/^ exp(— i/r^) 



(or OS harmonic oscillator wave function (f)ooo{r, u), which will be defined later) with 100 % 
occupancy (A = 1). This means that all particles are condensed in that single orbit, i.e. 
an ideal Bose- Einstein condensation is realized in the laboratory frame. This feature is 
independent of the number of particles N. 

Next we consider the internal one-particle density matrix for the A^-particle Bose-Einstein 
condensation in a harmonic trap described by the wave function Eq. (6) with the total Hamil- 
tonian Eq. (4) . Internal means that the density is free from the center-of-mass coordinate 
of the system, as it must be the case, when dealing with, e.g., selfbound systems. In the 
present paper, two kinds of internal coordinate sets are introduced; 1) coordinates with re- 
spect to the center of mass of the total system and 2) Jacobi coordinates. The former set 
was first considered by Pcthick and Pitaevskii [18] to define the internal onc-particlc density 
matrix. We call it Pethick-Pitaevskii-type (PP-type) internal one-particle density matrix in 
the present paper. For the latter set, we call it Jacobi-type density matrix. 



A. Pethick-Pitaevskii-type internal one-particle density matrix 



As already mentioned in the introduction, Pethick and Pitaevskii proposed to consider the 
internal single particle density matrix for a Bose condensed system when only the center-of- 
mass but no internal degree of freedom is excited. In order to define an internal one-particle 
density matrix, Pethick and Pitaevskii adopted internal coordinates defined with respect 
to the center of mass of the total system [18]. The center-of-mass coordinate R and the 
coordinate of particle i relative to the center of mass are given, respectively, by 



1 



N 



(10) 



Since J2f=i Qi = Oi only — 1 of the are independent. Hereafter, as the independent 
internal coordinates, we adopt q^ [i = 1, ■ ■ ■ , N — 1). We introduce the conjugate momenta 
TTi {i — 1, ■ ■ ■ , N — 1) and P for the coordinates q^ {i — 1, ■ ■ ■ , N — 1) and R, respectively. 
Then, the total Hamiltonian in Eq. (4) is rewritten as 



H = H, 



int 



rr 1 

" 2Nm 



N 



N-1 



i=l 



i<i'=l 



+ muj 



N-1 N-1 

E + E 9i- Qi' 

1=1 i<i'=l 



+ -Nmuj^R\ 



(11) 

,(12) 
(13) 



where ifint and Hem denote the internal and center-of-mass Hamiltonians, respectively. It 
should be mentioned that the above Hamiltonian includes the cross terms, • qr^/ and TTj • TTj/, 
which, in general, appear so far as one takes non-orthonormal coordinate systems [23-26] . 
By using the relation 

N N-l 



1=1 i,i'=l 

the total wave function in Eq. (6) is expressed as 

1 



*({n}f=i 



Ar3/2 



X $i,t({Q,}f=I^)*cm(i«), 



'N X (2i/) 



N-l\ 3/4 



TT 



N-l 



exp 



N-l 



iA'=l 



$em(i2) = i^-^) exp(-7VI/i^2), 



(14) 

(15) 
(16) 

(17) 



where $int and $cm denote the internal and center-of-mass wave functions, respectively 
The separability of $int and $cm comes from the fact that the total Hamiltonian is separable 
with respect to the internal and center-of-mass parts [see Eq. (11)] or directly from Eq. (14) 
inserted into Eq. (6). The factor l/N^/"^ in Eq. (15) drops out when taking into account 
the Jacobian d{{ri\f^]) / ^{{q^}^Sl^ , R) — N^, coming from the coordinate transformation 
from the space-fixed system to the center-of-mass system. It should be noted that both 
wave functions $int and $cm are normahzed and satisfy i^mt^int — (3/2) (A?" — l)^a;$int and 
-f^cm'^'cm = (3/2)/ia'$cm, respectively. 

According to Pethick and Pitaevskii [18], the internal one-particle density matrix is de- 
fined as, 

Pint,pp(9, q') = j dq^--- dq^_^ Pint,pp(g, ' " " , Qn-i: Q' ^ Q2r ■ ■ , Qn-i), (18) 
Pint,PP({g,}£l\ {q^l-i') ^N'JdR $*({q, + R}g,M{q', + R}g,), (19) 

where the Jacobian factor A^^ is inserted in the above equation, which is missing in Ref. [18]. 
Using the wave function $({?"i}ili) in Eq. (6), the internal one-particle density matrix (18) 
is expressed as 



N 



N-l 



3/2 /2z.\^/2 



TT 



exp 



3A^-2 , 2 ,2, N 

A<i + <i ) + 



2(A^- 1) 



N-l 



lyq ■ q 



(20) 



It is noted that this density matrix depends on the number of particles N and contains 
the cross term q • q'. The origin of the cross term comes from the nonseparability of the 



coordinate from the other {N — 2) coordinates {i — 2, - ■ ■ ,N — 1) in the internal wave 
function of Eq. (16), or equivalently from the fact that the internal Hamiltonian of Eq. (12) 
is non-separable between the first particle operators tti and q^ and the other (iV — 2) particle 
operators tTj and q^ {i = 2, ■ ■ ■ , N — 1). The result of Eq. (20) can also be verified with the 
use of the internal density matrix in Eq. (19) presented as follows: 



N-U 



(21) 



where ^*nt({9i}ilT^) is the internal wave function defined in Eq. (16). 

Let us discuss the nature of the internal one-particle density Pint(qf, q')- First we study 
the single-particle orbits and their eigenvalues obtained by solving the eigenvalue problem 
for the density matrix, 



j p^,,^AQ,Q')mdq'^M{q)■ 



{22) 



This equation can be solved analytically. There are several ways to solve it. In Appendix 
A we give one method. We can also analyze the description given in Ref. [27]. The single- 
particle orbits are expressed by the harmonic oscillator wave functions (f)nLM{Q, Pn) with 
the orbital angular momentum L, magnetic quantum number M and harmonic oscillator 
quanta Q = 2n + L (n = 0, 1, 2, • • •), 

"(2/5^)^+3/22"+^+2n!^'/' 



y/n{2n + 2L + l)\\ 



2A^ 



N 



! N-1 

The eigenvalues or occupation probabilities A are given as 



q'Li'^^'/'\2/3m')exp{-(3W)YLM{q), (23) 

(24) 



A, 



(LM) 
n,N 



AN 



3/2 



N -2 



-1 2n+L 



3N-2 + 2^2N{N-1)\ [sN -2 + 2^2N{N 
and satisfy the following completeness relation, 

oo L oo 

E E EAj^^^i. 

L=0 M=-Ln=0 



(25) 



(26) 



Then, the internal one-particle density matrix for the particles in Eq. (20) can be presented 
in terms of the wave functions (23) 



Pint,Pp(«'9') = E E 12 \(l>nLM{q:pN)) Xn,N^ {(l>nLM{q':pN)l 
L=0 M=-L n=0 



(27) 
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The occupation probability with respect to the partial wave with quantum number L is 
defined as 

L oo 

aS^^= E EA^r^. (28) 

M=-Ln=0 

In the macroscopic limit {N — > oo), the internal one-particle density matrix Pint,pp,Ar=oo' 
eigenfunctions 0int,pp,Af=oo) ^-nd eigenvalues \^n,N=oQ given by 



Pint,PP,Af=oo(9, Q) = [—) exp 



3/2 r 3 

-u{q^ + g") + i^q • q' 



(29) 

0mt,PP,7V=oo(g) = 4>nLM{q, V2u) , (30) 

A^^ioo = 2^(3 - 2V2)2"+^+3/2^ (31) 

L oo 

ASvioo= E EA^i2oo = (2^ + l)(2-V2)(3-2V2)^ (32) 

M=-L n=0 

We remark that the summed eigenvalues A^^^^^^ still depend on the angular momentum L. 

The eigenvalue or occupation probability of the PP-type internal one-particle density 
matrix, aS^^ of Eq. (28) depends on the orbital angular momentum L, and the particle 
number N. In the case oi N — 3, we obtain the occupation probabihties, A^'^^ = 0.804 for S 
wave, A^^^ = 0.173 for P, A^^^ = 0.021 for D, A^ ' = 0.002 for F, and so on. Increasing the 
particle number, the S'-wave occupation probability is decreasing, while the higher partial- 
wave ones are increasing (see Fig. 1 in Refs. [21, 22]). These results show that the PP-type 
one-particle density matrix leads to a fragmented condensate. The reason why the PP-type 
internal density matrix shows the fragmented condensate is due to the existence of the cross 
term q ■ q' in Eq. (20). 



B. Jacobi-type internal one-particle density matrix 

For the A^-particle system, we define the {N — 1) internal Jacobi coordinates {^j, {i = 
1, • • • , — 1)} and the center-of-mass coordinate R as follows: 

1 ^ 

^i = ^i-W^ ^ '''^^ for z = l,2,---,A^-l (33) 
^ k=i+i 

1 ^ 

1=1 
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where denotes the relative coordinate between the first particle and the remaining (N — l) 
particles, and other Jacobi coordinates are self-evident. Then, the A'"-particle Hamiltonian 
in Eq. (4) can be separated into the internal and center-of-mass Hamiltonian, 

H — i?int + -f^cm, (35) 

^-1 1 1 N-i 

i/em = -^P^ + -Nmw^R^, (37) 
^ 2Nm 2 ' ^ ' 

where TTj and P denote the conjugate momenta corresponding to the coordinates and R, 

respectively. 

Since the total Hamiltonian is a sum of decoupled N harmonic oscillator Hamiltonians 
( for R and (A^ — 1) internal Jacobi coordinates ), the total wave function in Eq. (6) is 
expressed as 

H{ri}l,) = ^int({£.}£l^)$cm(H), (38) 

/o,, X 3/4 

$.nt({uiii^) = n (v) ^-pH^?) (39) 

i=l ^ " ^ 



$em(i?) = [—^) expi-NiyR'). (41) 
This expression of ^{{ri}^i) can also be derived by inserting the relation 

Erl^NR^+Y:'^^ (42) 



=1 i=i ^ 



into the original expression given in Eq. (6). The total wave function is the product of the 
internal and center-of-mass wave functions, $int and $cm, which are normalized and satisfy 
the relations, -f^mt'^'int — (3/2)(A'" — l)^a;$int and -ffcm^cm — (3/2)^<l>cin- It is noted that the 
internal wave function is given as a product of harmonic oscillator wave functions 0ooo(£i, J^i) 
(i = l,---,7V-l). 

The Jacobi-type internal density matrix writes 
Then, the Jacobi-type one-particle density matrix is defined with respect to and ^'i as 

Pinij(^> = / • • • d^N-l Pint,j(^, ^2, ' ' ' , ^N-l, ^' , ^2, ' ' ' , ^N-l), (44) 

12 



N 



N 



N 



(45) 



This choice of the coordinate for the internal density matrix is natural, because the 
single particle orbit should be defined with respect to the relative coordinate between one 
particle and the other remaining — 1 particles in the Jacobi coordinate system. The 
compatibility between the bosonic symmetry of the system and the above definition of 
the internal one-particle density matrix, where it may seem that one special coordinate is 
singled out, is explained in Appendix B. The eigenvalue equation of the one-particle density 
matrix can easily be solved analytically. We find that the density matrix has only one non- 
zero eigenvalue A = 1 with corresponding eigenfunction 0int,j,Af which is the 05' harmonic 
oscillator wave function with 100 % occupancy (A = 1), 



0int,J,iv(C) = (^LO(^MO0OOO(€, (A^ - l)i^/A^), 

>(LAf) 
At 



(46) 
(47) 

(48) 



M=-L 

This means that all particles are condensed in the single 05 particle state, although the 
size parameter in the state [Eq. (46)] depends on A^ and is slightly different from that in 
the eigenfunction 0ooo(^) ^) in laboratory frame, discussed in the beginning of Sec. II. In 



the macroscopic limit (A^ oo), the internal one-particle density matrix p\J^ 



(1) 



,PP,iV=oo 



eigenfunction 0int,pp,Af=oo, and eigenvalues \^n,N=oo given by 



Pmt,J,N= 



(tO = (^) exp 



TT 



0int,J,iV=oo(€) — ^LqSmq(PoOo{$,^) 



(LM) 
N=c 

M=-L 



A 



JLO- 



, its 

(49) 
(50) 
(51) 
(52) 



The Jacobi-typc onc-particlc density matrix for finite particle number is separable with 
respect to the coordinates ^ and This leads to the only one non-zero eigenvalue A = 
1 which is the same as for an ideal Bose-Einstein condensate system in the laboratory 
system. Also the internal density matrix (49) becomes identical to the one in the laboratory 
system, Eq. (8), for N — oo. The separability of the internal density matrix comes from 
the separability of the internal wave function in Eq. (40) with respect to different Jacobi 
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coordinates, which originates from the separabihty of the internal Hamiltonian in Eq. (36) 
with respect to different Jacobi coordinates. This feature has it origin in the fact that the 
Jacobi coordinates form an orthogonal coordinate system. 

III. CRITERION FOR THE CHOICE OF ADEQUATE INTERNAL COORDI- 
NATES 

A. Convergence to Bose-Einstein Condensation in the Macroscopic Limit 

In the previous section we learned that the outcome of the diagonalization of the internal 
density matrix depends on the choice of the internal coordinates. This is a serious problem 
for treating condensates in internal self-bound systems such as a-particle condensates in 
nuclear systems or small droplets of superfiuid ^He, because only internal degrees of freedom 
are relevant in these systems. 

In order to overcome the difficulty, we gave a criterion for the choice of the internal 
coordinates: In the macroscopic limit {N — > oo) the internal density matrix should have 
the same eigenvalues and eigenf unctions as those of the ideal Bose-Einstein condensate in 
the laboratory frame. This is a very physical boundary condition. For understanding this 
physical condition, consideration of the following situation in the macroscopic limit may 
be helpful: In the laboratory frame, the center-of-mass motion of the present system is 
described by the wave function $cm(-R) of Eq. (17). In the macroscopic limit, the center- 
of-mass coordinate should be at the coordinate origin of the laboratory frame, because the 
probability of finding the center-of-mass coordinate at position R is given by |$cm(-R)P = 
(2A^i//7r)^/^ exp(— 2A^i/il^) which becomes a delta function S{R) in the limit of — > oo. 
Thus the internal coordinate should have the same meaning as the position of the i-th 
particle coordinate in the laboratory frame. Also the first Jacobi internal coordinate 
should have the same meaning as the position of the first particle coordinate ri because 

= {N/{N - l))qi ^ n for AT ^ oo. 

As shown in Eqs. (30) and (31), the PP-type internal density matrix in the macroscopic 
limit does not satisfy the condition, while the density matrix of the Jacobi-type fulfills this 
condition as given in Eqs. (50) and (51). As already mentioned, the reason why the PP-type 
internal density matrix does not satisfy the physical boundary condition is due to the fact 
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that it exhibits a nonvanishing cross term or correlation term q ■ q' in the case oi N — oo, 
originating from the pseudo two-body interaction terms in the internal Hamiltonian. On 
the other hand, in the Jacobi-type internal density matrix, the pseudo two-body interaction 

terms disappear. 

These results mean that one should take internal coordinates which do not produce any 
correlation in the internal one-particle density matrix in the macroscopic limit. Otherwise, 
unphysical situations occur like for the PP-type internal density matrix. One choice fulfilling 
the physical condition is the one of the internal Jacobi coordinates. Of course, there are 
many sets of internal coordinates which satisfy the physical condition. For example, they 
may be those satisfying the following two conditions: (1) The relative coordinate ^ between 
one particle and the other remaining particles should be used as the coordinate of the internal 
one-particle density matrix p^^^, and (2) the internal one-particle density matrix be separable 
with respect to ^ and ^' in the macroscopic limit, p-nt (2i//7r)^/^ exp[— for 

N ^ oo. Among the coordinates satisfying the two conditions, Jacobi coordinates are very 
convenient and useful to describe the internal Hamiltonian and thus the internal one-particle 
density matrix. We will see in the next section that this is in fact the only choice. 

B. Orthogonality of the Coordinate with respect to the {N — 2) Internal Coor- 
dinates 

In the previous section, a criterion was discussed for the choice of the adequate coordinate 
systems of the internal density matrix. Here we present the other criterion that the coor- 
dinate of the internal one-body density matrix, should be orthogonal to the remaining 
(A^ — 2) internal coordinates, where the (A^ — 2) coordinates do not need to be mutually 
orthogonal. This is called "maximal independence of the coordinate from the rest of 
(A" — 2) internal coordinates" . 

The coordinate of the internal one-particle density matrix should be either the one particle 
coordinate measured from the total center-of-mass or the one particle coordinate measured 
from the center-of-mass of the other (A^ — 1) particles. The PP-type internal one-particle 
density matrix adopts the former type coordinate qi while the Jacobi-type internal one- 
particle density matrix adopts the latter type coordinate ^i, coinciding in the macroscopic 
number. These two types of coordinates are essentially the same since they are related by 
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Qi — ii-^ ~ constructing the internal one-particle density matrix, we inte- 

grate the internal density matrix over the other {N — 2) internal coordinates. In order for 
the internal one-particle density matrix to have maximum information on the one-particle 
dcgrce-of-freedom of the system, these {N — 2) internal coordinates should be maximally 
independent from the coordinate or Since this requirement seems natural, we will 
study in this subsection its consequence. We will see below that this requirement implies 
the criterion proposed in the preceding subsection. 

The requirement that the coordinate or is maximally independent from the other 
(N — 2) internal coordinates means mathematically that coordinate qi or should be 
orthogonal to the other {N — 2) internal coordinates. The Jacobi coordinates are just such 
coordinates, since they constitute an orthogonal coordinate system. On the other hand 
the internal coordinates {qfj, i = 1 ~ (A^ — 1)} do not satisfy this requirement since they 
constitute a non-orthogonal coordinate system. For the sake of self-containedness we recall 
here the meaning of orthogonality between coordinates. Two coordinates /3 and 7 are 
defined to be mutually orthogonal when their expansion coefficients of linear combination 
with respect to the A'" particle coordinates {r^, (i = 1 ~ N)}, {Ci(/3), (i = 1 ~ N)} and 
{Q(7), (z = l-iV)}, 

N N 

f3 = J2Qif3)n, l = J2C^h)r^, (53) 
i=i 1=1 

are mutually orthogonal, J2iLiCi{/3)Ci{^) = 0. This definition can be stated as fol- 
lows. To any coordinate d we associate an A^-dimensional number vector C{5) like 
C(/3) = {cm, = 1 - N)} for f3 and C{j) = {^(7), = 1 ~ N)} for 7. If (3 7, 
C{^) 7^ C(7). ^ and 7 are said to be orthogonal when C(/3) • C(7) = 0. Needless to say, 
the total center-of-mass coordinate R is orthogonal to any kind of internal coordinates ^int, 
i.e. C{R) ■ C(5int) = 0. 

Let an internal coordinate system {ry^, {i = 1 ^ {N — 1))} be such a system satisfying 
the above-mentioned requirement. The coordinate rj^ is either q^ or and hence we here 
fix it as = $1 . The orthogonality of to the other (A^ — 2) coordinates {77^, {i = 
2 ^ {N - 1))} means C(|i) • C(77j = 0, {i ^ 2 ^ {N - 1)). Since the internal Jacobi 
coordinate system {i — 1 (N — 1))} is an orthogonal coordinate system, there the 
relations C{^i) ■ C(^j) = 0, {i — 2 {N — 1)) holds. Therefore the subspace spanned by 
the A'"-dimensional vectors {C(r7j), (i = 2 ~ (A?" — 1))} is identical to the subspace spanned 

16 



by A'"-dimensional vectors {C(^j), {i — 2 {N — 1))}. We, therefore, have 



N-l 



2~ (A^-1)). 



(54) 



This relation is equivalent to the relation 



N-l 



{N-l)). 



(55) 



By inserting Eq. (55) into Eq. (42) we obtain 



N 




(56) 



1=1 




(57) 



i,i'=2 




E dj,idj,i'- 
j=2 



(58) 



It is to be noticed that unlike the Jacobi coordinates the mutual orthogonality within the 
(A^ — 2) coordinates {rj^, i = 2 ~ (A^ — 1)} does not hold in general. 

By using Eq. (57) we can calculate the internal one-particle density matrix as follows 



This result shows that Pint ('?)'?') is just the same as the Jacobi-type internal one-particle 
density matrix. Thus we see that our above requirement concerning the internal coordinates, 
gives us a unique result for the internal one-particle density matrix. It is just identical to 
the Jacobi-type internal one-particle density matrix. 

In the previous section III A we have proposed also another requirement, namely that the 
internal one-particle density matrix should converge to the one-particle density matrix in the 
laboratory frame in the macroscopic limit. As we discussed in Sec. Ill A, this requirement 
implies that there should not appear any cross terms of with the other internal coordinates 
in the internal wave function at least in the macroscopic limit. Clearly the absence of the 
cross terms of with the other internal coordinates is realized only when the coordinate 




(59) 



(60) 
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£i is orthogonal to all the other {N — 2) internal coordinates. Thus the criterion in the 
preceding subsection results from the present requirement of the orthogonality of to the 
other (A^ — 2) internal coordinates, where the {N — 2) coordinates generally do not need to 
be mutually orthogonal. 

One may argue that a physical quantity should not depend on the choice of the coordinate 
system. In the case of the one-particle density matrix this argument can be true under the 
condition that we have extracted maximum information on the one-particle degree of freedom 
of the system. As already noticed this condition is the same as the requirement that the 
coordinate of the one-particle density matrix is orthogonal to all the other coordinates of the 
system. This is also true in the macroscopic system where wc usually adopt the coordinate 
system composed of individual particle coordinates (rj, (i = 1 ~ N)}. The coordinate Vi 
of the one-particle density matrix is of course orthogonal to all the other coordinates. If 
we adopt the coordinate system {R, q^, [i — 1 [N — 1))} and calculate the one-particle 
density matrix with respect to the coordinate which is practically the same as fi in the 
macroscopic system by integrating out the remaining coordinates {R, q^, (i — 2 (N — l))}, 
we get the result given in Eq. (29), namely (2z//7r)^/^ exp[—{3/2)u{ql + q'^^) + uqi ■ q[]. This 
inadequate result is, of course, due to the non-orthogonality of qi to the other coordinates 
{qfj, {i — 2 {N — 1))}, as we have already seen earlier. 

The PP-type internal one-particle density matrix Pint,pp different from the Jacobi-type 
one Pint^j- This is true for the non-diagonal elements but for the diagonal elements they are 
the same except for the Jacobian factor d{qi) / d{^i) — {{N — 1)/A^)^, (which goes to unity 
in the macroscopic limit) 

/ N 

Piipp(9i'^i) = (61) 
It is to be noticed that both Pint,pp a-nd p{^^ j are normalized, 

/ dq.pSlMli^ 9i) = / dLpSlA^i: L) = 1- (62) 

The equality of Eq. (61) is easily proved by using Eqs. (20) and (45). Let us study these 
relations between p^^]. pp and j a little more in detail. First we note the relation between 
the two coordinate systems, {qfj, (i = 1 ~ (A^ — 1))} and {^j, {i — 1 {N — 1))}, 

1 N -1 

= - g ivTw + NTI^fi- ^^^^ 
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When we calculate pl^^ j, we make the following integration 

^'l) = / ^^^2 • • • ^^^AT-l ^int(^l, £2, • • • , ^N-l)^U^'l: ^2, " " " , ^N-l) (64) 

x^mt{q[,q2,---,Q'N-i) (65) 

= 9.(^1 -ri) = g. + ^(li-ri), (66) 

where 

^int(qi, qf2, • • • , Qn-i) = ^int(^l, £2, ' ' ' , £iV-l), (67) 

^int(q'l, 92, • • • , q^-l) = *mt(^'l, ^2, ■ ■ ■ , ^N-l)- (68) 

For non-diagonal elements of Pint,j(£i, ^'i), since 7^ it follows that q- 7^ q^ for all 
i — 1 (N — 1). This tells us that Pint,j(^i, £'i) can not be proportional to Pint,pp{Qii Qi) 
because the latter is obtained from 



/ 



dq2 • • • dq^_i ^int(qi, 92, • • • , qN-i) ^int(qi, ^[2, ' ' " , 9jv-i)- (69) 



On the other hand, for diagonal elements, Pint,j(£i, ^i), since = ^'i, there follows q- = q^ 
for all i = 1 ~ (A?" — 1). This tells us that Pint,j(^i, ^i) is now proportional to Pint,pp(Qi, 9i)- 
The appearance of the term (1/A^)(^^ — $,[) in the relation q'- = q, + (l/iV)(^^ — ^[) just 
stems from the non-orthogonality of qi to q^ (i = 2 ~ (A^ — 1)). 

We have proved in this section that the internal one-body density matrix is uniquely 
determined for the OS harmonic oscillator wave function. However, it is noted that this 
uniqueness holds in the case of a general wave function including the 05' harmonic oscillator 
one. In fact, Suzuki et al. [29, 30] have already given the proof of the uniqueness for the 
more general wave function ^ which is expanded in terms of the correlated Gaussian basis 
g. The explicit forms of ^' and g are expressed as follows: 

^^Y.CkSg{s^'^;A^'\x) (70) 

k 

g = g{s; >1, a;) = exp [ - ^ ^ AijXi ■Xj + Y.Si- , (71) 

id i 

where x — {aji, a;2, • • • , Xn-i} and S are a set of internal coordinate of A'"-boson system and 
the symmetrization operator acting on the N bosons, respectively. C^, A'^^'' and s^^'> are the 
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expansion parameters, in which A^^^ denotes a symmetric positive-defined (N — l)x{N — 1) 
matrix, and s^'^^ represents s^'^^ — {si'\ ■ ■ ■ , s^^'Li} . The correlated Gaussian basis g 
is often used in ab-initio calculations and has succeeded in describing structures of many 
few-body systems [31]. The uniqueness of the internal one-body density matrix p for the 
wave function ^ in Eq. (70) is proved [29, 30] under the condition that one takes the 
following set of coordinates y — {l/i, 2/2; ' " " ? 2/iv-i} (obtained by a linear transformation 
from the coordinates x) as the coordinates of ^ adopted in calculating p: is chosen as the 
coordinate (t/i) of the internal one-body density matrix and the remaining {N — 2) internal 
coordinates (t/2, 2/3, • " " > Vn-i) orthogonal to q^, although the (iV — 2) coordinates do not 
need to be mutually orthogonal. This requirement for the coordinates y just corresponds 
to the "maximally independence of the coordinate (q^) from the rest of {N — 2) internal 
coordinates" as mentioned above. The proof we gave in this section is specified to the simple 
system which is composed of N bosons in the harmonic oscillator potential and has been 
discussed by many authors (for example, see Refs. [18, 19, 21, 22]). 



IV. SUMMARY 



We investigated the internal one-particle density matrix in the case of ideal Bose-Einstein 
condensates with a finite number ( iV ) of particles in a harmonic trap. Wc calculated the 
explicit form of the internal one-particle density matrix following the definition of Pethick 
and Pitaevskii (PP) and solved its eigenvalue problem. The result was found to show a 
fragmented condensate, contrary to what PP expected. On the other hand the Jacobi- 
type internal one-particle density matrix gives us complete condensation. It means that the 
internal one-particle density matrix is different in general for different choices of the internal 
coordinate system. In this paper we outlined two physically motivated criteria for the choice 
of the adequate coordinate system leading to a unique answer for the internal one-particle 
density matrix. One criterion is that in the infinite particle number ( N — 00 ) limit the 
internal one-particle density matrix should have the same eigenvalues and eigenfunctions 
as those of the corresponding ideal Bose-Einstein condensate in the laboratory frame. The 
other criterion is that the coordinate of the internal one-particle density matrix which is 
either q^ = Vi — R or = {N/{N — 1))^!, should be maximally independent from the 
remaining [N — 2) internal coordinates. Mathematically this criterion means that (or £^) 
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is orthogonal to the remaining {N — 2) internal coordinates, though the {N — 2) coordinates, 
in general, do not need to be mutually orthogonal. This second criterion was shown to 
imply the first criterion. We saw that the internal Jacobi coordinate system satisfies these 
two criteria while the internal coordinate system adopted by Pethick and Pitacvskii for the 
construction of the internal one-particle density matrix does not. Furthermore we argued 
that these two criteria uniquely determine the internal one-particle density matrix which is 
the same as that calculated with the Jacobi coordinates. The results of this paper justify 
the use of the Jacobi coordinates in Refs. [16, 17] where the Bose condensation of a few 
a-particles in extended states of self-conjugate light nuclei was considered. However, our 
results arc of more general interest. For example the number of bosons captured in each 
site of an optical lattice is often very small [28] and, therefore, our analysis surely applies 
to that situation as well. We also believe that our present study has considerably clarified 
the somewhat controversial issue of how to define a Bose condensate of a finite number of 
particles in their internal coordinate system. This question is particularly, but not only, 
relevant for selfbound Bose systems, like a loosely bound gas of a-particles in nuclei or 
nano-droplets of liquid ^He. 

The authors thank Kanto Gakuin University (KGU) and Yukawa Institute for Theoretical 
Physics at Kyoto University, Japan. Discussions during the KGU Yokohama Autumn School 
of Nuclear Physics and the YIPQS international molecule workshop held in October 2008 
were useful to complete this work. 

APPENDIX A: SOLUTION OF THE EIGENVALUE PROBLEM OF THE SYM- 
METRIC GAUSSIAN INTEGRAL KERNEL 

The one-particle density matrices discussed in this paper are symmetric Gaussian integral 
kernels of the form 



{r\Q\r')^Q{r,r') 




exp{— a(r^ -|- r'^) -|- br ■ r'}. 



(Al) 



This kernel is normalized 




(A2) 



for which we need the condition 2a — 6 > 0. 
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The eigenvalue problem of the kernel Q{r,r') can be solved analytically. There are several 
ways to solve it. We can e.g. use the method given in Ref. [27]. Here we explain another 
procedure which gives us the following answer (see below) 

I/O/ w ocloc. 



Here al and dc are creation and annihilation operators of harmonic oscillation of size pa- 
rameter c, respectively. The eigenfunctions of Q are just the same as those of the operator 
ocl ■ cXc which are of course the harmonic oscillator functions of the size parameter c. Thus 
we have 

Q|0nLM(c)) = Ai^'^^)|0„iM(c)), (A5) 

{r\(f)„LMic)) = (pnLM{r, c) , (A6) 

The explicit form of (f>nLM{r,c) is given in Eq. (23). It is to be noted that the eigenvalue 
^(L,M) Jepgnds only on the number of harmonic oscillator quanta (2n + L). It is just the 
SU(3)-scalar property of the operator Q. A merit of our method as explained below, is that 
it gives us a clearer understanding of the SU(3) symmetry of the operator Q. 

The relation of Eq. (A3) comes directly from the following theorem [32] valid for a general 
operator P: 

P^Kp{p{ci\,a^)}, (A8) 
where Afop is the operator of normal ordering, and 



«7 



= exp{-Z* ■ Z'){A^{Z)\P\A^{Z')), (All) 



{r\A,{Z)) = A,{r,Z) = exp 



ZV 1 



(A12) 
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The state \A-y{Z)) is the well-known coherent state of harmonic oscillator of size parameter 
7, 

\A^{Z)) = exp{Z-ay \A^{Z = 0)), 
a^\Aj{Z)) = Z\A^{Z)), 
A^{r,Z) = A^{x,Z^)A^{y,Zy)A^{z,Z,), 



Aj[x, Zx) — 



27\4 
— exp 

TT / 



(A13) 
(A14) 

(A15) 

(A16) 



n=o vn! 



[a. 



t )n 

■y,xJ 



-Aj[x, Z^ — 0). 



(A17) 
(A18) 



The proof of the general theorem of Eq. (A8) is quite easy. It is given by comparing the 
matrix elements of two operators, P and N'op{p{oct^, oc^)}, formed with the coherent states 
\Aj{Z)) which constitute an overcomplete set of states. This is shown by the following 
operation: 



{A,{Z)\KMccl,cc,)}\A,{Z')) = p{Z\Z'){A,{Z)\A,{Z')) 

= {A,{Z)\P\A,{Z')). 



(A19) 
(A20) 



The equahty of Eq. (A19) is due to Eq. (A14). 

By applying the above general theorem to our present operator Q of Eq. (Al), we obtain 

{A,{Z)\Q\A,{Z')) 



q{Z\Z) 



{A,{Z)\A,{Z')) 

^2 



(7 + a)2-(62/4), 



exp 



F(Z*^ + Z'^) + GZ* ■ Z' 



1 



F = -- + 
G = -1 + 



7(7 + a) 



2 (7 + a)2 - (62/4) ' 
76 

(7 + a)2 - (62/4) ■ 



By choosing 



^ — c — \ a'^ 



62 

T' 



we obtain F = 0. Thus we have for 7 = c 

,(Z-,Z')=f^)*exp 

\ a + c 



2(a-c] 



1\Z* ■ Z' 



(A21) 

(A22) 
(A23) 
(A24) 

(A25) 

(A26) 
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By using the well-known formula [33] 

M>p{exp[(r - l)at . a,]} = (A27) 
we obtain the desired result of Eq. (A3) 

Q = Kp{q{(xl,(x,)} (A28) 

The formula of Eq. (A27) can also easily be proved by comparing the matrix elements of 
A/'op{exp[(y — l)Q;jy • ocj]} and Y^^t''-^'' formed with coherent states. For treating Y^^-t'^^, 
the expansion of the coherent state by harmonic oscillator functions given in Eq. (A17) is 
useful. 



APPENDIX B: BOSONIC SYMMETRY OF THE JACOBI-TYPE INTERNAL 
ONE-PARTICLE DENSITY MATRIX 

The bosonic symmetry of the Jacobi-type internal one-particle density matrix defined in 
Eq. (44) can be seen in the following identity relation, 

O = / • • • d^^_, $int(e, ^2, • • • , ^iV-l)^int(r, ^2, " " " , ^AT-l) (Bl) 

= j d^id^2 ■ ■ ■ d^N-i X ■ ■ ■ d^N-i 

X^int (£1, ^2, • • • , ^N-l)0{l, t 0*int(£'l, ^2, ' ' ' , C-l) (B2) 

= j d^id^2 ■ ■ ■ d^N-i X d^[d^2 ■ ■ ■ 

1 N 

X$int(£l, ^2, • • • , £iV-l)T^ E 0{k, C 0^int(^'l, ^2, ' ' ' , ^N-l), (B3) 

fc=l 

o(fc,to = 5(€i-€)^(€f-on^(^j-^f)- (B4) 

Here we define N different sets of Jacobi coordinates (^i, ^2:' " ■> ^n-i) by N cyclic permu- 
tations of particle indices (1, 2, • • • , A?") 

1 ^ 

= rpk{i) - T7— E ^Pk{jh (B5) 
{pk{l),Pk{2), ■ ■ ■,PkiN)) = {k,k + l,---,N,l,2,---,k-l), {k = lr.N). (B6) 
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We should note that, since $int(^i,£2) ' ' ' j^n-i) totally symmetric for particle permu- 
tations as is clear from the relation $(ri, r2, • • • , Vn) — $int(^i, ' " ■> ^Ar-i)^cm(-R), we 
have 

^int(£l, ^2, • • • , ^AT-l) = ^Utl, £2, • • • , ^^-1) (k^lr^N). (B7) 

And of course we have 

d^,d^2 ■ ■ ■ c?£iv-i = de,de2 ■ ■ ■ deN-,. (B8) 

The expression of Eq. (B3) is the manifestation of the bosonic symmetry of the Jacobi-type 
internal one-particle density matrix. This equality is due to the following relations 

{0{k, 1 10) = (0(1, t O), {k = 2^N),, (B9) 

where {0{k, ^, ^')) is defined as 

{0{k, ^, ^')> = / d^,d^, ■ ■ ■ d^j,_, X d^[d^', ■ ■ ■ d^'j,_, 

X^int (€1, I2, • • • , ^N-l)0{k, t €')*int(€'l, ■ ■ ■ , C-l), (BIO) 

(A; = 1 ~ N). 

We can prove the equality of Eq. (B9), by using Eq. (B7) and Eq. (B8), as follows 
(0(fc,tO) = Jd$,d^,---d^j,^,xdi[d^',---d^'j,_, 

X$int(£l, ^2, • • • , ^N-l)0{k, I, O^intd'l, £2, " " " , £'iV-l) (BH) 

= j d^id^2 ■ ■ ■ d^%-i X d^'l'd^2 ■ ■ ■ d^N-i 

x<f i„t (^t ■ ■ ■ , eN-i)m t O'^iA^i, €2', • • • , €^-1') (B12) 

= / c?€ic?€2 • • • d^N-1 X d^[d^2 ■ ■ ■ d^'^_i 

X$int(^l, ^2, • • • , ^iV-l)<5(l, t 0^int(^'l, ^2, ' ' ' , ^^-l) (B13) 

= {0{l,ta- (B14) 
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